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An update is given of the determination of leptonic decay constants of charm 
and beauty mesons in the framework of relativistic Hilbert moments and Laplace 
transform QCD sum rules. 
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In this talk I present an update of relativistic QCD sum rule estimates of the leptonic 
decay constants 

(0\A^\P(k))^iV2f P k^, (1) 

(P = D, D s , B, B s ), so that in this convention f n = 93.2 MeV. The current theoretical 
status of this problem is rather confusing, due to the existence of a plethora of predictions 
not always in mutual agreement, even if obtained in the same framework. The major 
difference between the various determinations to be found in the literature |l|, ^ |3|, f|, || is 
due to the particular choice of input parameters in the sum rules. Among these, the values 
of the heavy quark mass and the asymptotic freedom threshold impact the most on the 
prediction for fp. To a lesser extent, sum rule windows or stability criteria have also an 
influence on the results. Hence, it is very important to keep this in mind when comparing 
different determinations of the same quantity. It is then not very illuminating to present 
a table showing all the various existing predictions. Instead, I shall present a single 
prediction for each fp, from Hilbert moments and, separately, from Laplace transform 
QCD sum rules, using the criterion that the experimental value of the meson mass should 
be reproduced by the sum rules, when implemented by our current best knowledge of the 
input parameters. Considerable progress has now been made in improving the accuracy 
and reliability of the determinations of quark masses |§, [[?]], and vacuum condensates 
|| entering the QCD sum rules. The uncertainties in these parameters reflect in the 
uncertainty on fp, for a given type of sum rule. The two types of sum rules do not give 
exactly the same answer for a given fp, although the results are not all that different. 
A comparison between the two separate determinations can provide a feeling for the 
systematic uncertainties involved in this approach. 
In order to estimate fp one considers the two-point function 

ij 5 (q 2 ) =i j d A x e^(0|T (&> A^x) d v 4(0)) |0> , (2) 

where <9 M A^(x) = (ttiq + m q ) : q(x)ij 5 Q(x): with q{Q) being the light (heavy) quark 
field and m,(mg) its corresponding QCD (current) mass. The function ^(Q 2 ); Q 2 — ~<? 2 ; 
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satisfies a dispersion relation 
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<Is — ^LJ- + subtractions, 



(3) 



s + Q 2 

defined up to two subtractions, which can be disposed of by taking at least two derivatives 
in (3). In this fashion one obtains the Hilbert power moments, which at Q 2 = become 



[n 
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The point Q 2 = is appropriate for heavy-light quark currents, to the extent that (p n can 
be computed in perturbative QCD, adding non-perturbative power corrections which fall 
off by powers of the heavy quark mass. These corrections are parametrized by vacuum 
expectation values of the quark and gluon fields in the QCD Lagrangian, and are organized 
according to their dimension. For instance, in the limit m q — > 0, well justified for D u & 
and B u d mesons, the perturbative contribution to (p n (0) to order 0(a s ) is given by 

o / , \n-l 

V (n) (0)\ PT = -U) B(n, 3)[1 + a° n a s ], (5) 



where B(x,y) is the beta function, a s 
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The non-perturbative part, always in the limit m q — > 0, becomes 
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where p is a measure of the deviation from the vacuum saturation approximation of the 
four-quark condensate (p\vs — !)■ I n the case of the D s and S s mesons, where the 
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approximation m q = should not be made, the full expressions given in |J must be 
used for Lp^(0). Finally, the hadronic spectral function appearing on the r.h.s. of (4) 
is parametrized by the ground state pseudoscalar meson pole plus a continuum starting 
at some threshold Sq. This continuum is expected to be well approximated by the QCD 
spectral function, computed in perturbation theory, provided So is high enough, i.e. 

-lmi/) 5 (s)\HAD = VpMf, 5{s - m%) + 9{s - s )- lm^ 5 {s)\ PT • (8) 

7T 7T 

By taking the ratio of any two consecutive moments one obtains an expression for Mp as 
a function of sq, the latter being a - priori unknown. The calculation will be meaningful 
provided M P does not depend strongly on s , i.e. there should be a relatively wide range 
of values of So leading to a value of Mp with a reasonably small spread. This is certainly 
the case for D and D s , where one obtains using the first two moments (n = 1, 2) 

s = 2M%-3M% , M D = 1.85 ± O.lSGeV , (9) 

s = 2M 2 Ds - 3M 2 Ds , M Ds = 1.9 ± O.lGeV , (10) 

to be compared with the experimental values: M D \ EX p = 1-87 GeV, and M Ds \ EX p = 1.97 
GeV. With increasing heavy quark mass, the perturbative contribution increases in impor- 
tance relative to the non-perturbative part. Therefore, the stability region in s becomes 
narrower. For instance, for Q=b one finds that with sq ~ (1.1 — 1.3)M B , the predicted 
mass is Mb = 5.2 ± 0.2 GeV (Mb\exp = 5.27 GeV). An additional confidence criterion 
often invoked is that of the hierachy of the non-perturbative power corrections. An in- 
spection of Eq.(7) shows that this hierarchy is not respected in the case of Q=c, as the 
dimension d=5 term can easily become bigger than the d=4 term. In addition, the d=5 
contribution could become bigger than the perturbative contribution. In fact, there are 
no real solutions for f D and fc a , unless Mq < 0.5GeV 2 , and n < 2, as noticed in 0. 
These problems are absent for Q = b, where the perturbative term dominates the sum 
rule. In this case there are real solutions for all values of n, and M5 2 could be as high 
as 1 GeV 2 . In summary, there are certain advantages and shortcomings of the Hilbert 
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moments, which should be kept in mind when comparing results in this framework with 
those from Laplace transform sum rules. 

The criterion I shall adopt here is to fix so in such a way as to reproduce the experimental 
value of the pseudoscalar meson mass, for a given set of input parameters. The latter 
are chosen as follows: m c = 1.35 ± 0.05 GeV, m b =4.72 ± 0.05 GeV, m s ~ 190 MeV, 
A = 200 ± 100 MeV, p = 3 ± 1, < a s G 2 >= 0.038 - 0.11 Gel/ 4 , < qq > \ c = -0.010 
Gel/ 3 , < qq > \ b = -0.014 GeV 3 , and M 2 = 0.5 - 1.0 GeV 2 , except for Q=c where 
Mq = 0.5 GeV 2 , as mentioned above. Concerning f Dg and /b s , previous analyses |^|, |4[] 
have been made keeping m s ^ in the perturbative part, but not in the non - perturba- 
tive expression of the two-point function (2). This can be improved by keeping m s ^ 



everywhere (details are given in ||10|| ). An additional improvement, in the case of fs B is 



possible thanks to the recent measurement of the B s mass ||1 1|| : Mb s \exp — 5.37 GeV. 
The results from these Hilbert moment QCD sum rules are shown in Tables 1 and 2, 
for the two extreme values of the corresponding heavy quark masses. The minimum and 
maximum values of the leptonic decay constants are obtained by varying the input pa- 
rameters within the limits given above. Of particular importance are the results for the 
ratios /d,//d 5 and fsj Ib, which are basically independent of the heavy quark mass, and 
quite stable against changes in the rest of the input parameters. 



Table 1: Hilbert moment QCD sum rule results for Q=c 





m c = 


1.3 GeV 


m c = 


1.4 GeV 




1.29 


- 1.79 


1.06 


- 1.51 


s \d (GeV 2 ) 


4.5 


- 8.0 


4.0 


- 6.0 


fD a 1 fir 


1.59 


- 2.06 


1.33 


- 1.74 


s \d 3 (GeV 2 ) 


5.0 


- 7.5 


4.5 


- 7.0 


IdJ Jd 


1.23 


- 1.15 


1.25 


- 1.15 



It is possible to choose a different kernel in the dispersion relation (3) and obtain other 
types of QCD sum rules, e.g. the Laplace sum rules 
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Table 2: Hilbert moment QCD sum rule results for Q = b 





m b = 4.67 GeV 


m b = 4.77 GeV 


SbI fir 


1.19 - 1.40 


1.02 - 1.18 


s \b (GeV 2 ) 


32.5 - 34.0 


32.0 - 33.0 


/b s / fn 


1.48 - 1.68 


1.26 - 1.44 


s \b s (GeV 2 ) 


34.0 - 35.0 


33.0 - 34.0 


SbsI fs 


1.24 - 1.20 


1.24 - 1.22 



2f 2 M P exp(-Mp/M 2 ) = / ds exp(-s/M 2 )-Iimp 5 (s)\ QC D 



+m 2 Qexp(—m 2 Q / M 2 ) 



< (a s /12n)G 2 - m Q qq > -±^§ (l - 2M 2 < qq > 



-i / 1 4 

1 / m n m 



2 -77TH-7n7Z (16/9)vr« s p < qq >< 



6M 2 V 2M 2 6M 4 
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In (11) m s = is understood (for details of the case m s ^ see ||10|). Notice that the sign 
of the gluon condensate is positive, contrary to that in [|IJ which is incorrect. Taking the 
first derivative with respect to M 2 in (11) gives an independent sum rule which can be 
used together with (11) in order to get an expression for the meson mass Mp, independent 
of fp. This procedure, which essentially fixes s , is quite important, i.e. a determination 
of fp will be reliable provided that Mp comes out right. This point has not been fully 
appreciated in some of the existing analyses. 

It has often been claimed that Laplace transform QCD sum rules are superior to the 
Hilbert moments for the determination of fp. I do not quite agree with this claim. 
In fact, in spite of the exponential kernel in the dispersion relation, results are rather 
sensitive to so- For instance, in the case of Q=c, a 20% change in s , around the value 
which gives the correct meson mass, induces typically a 10 % change in that mass as 
well as in fp,. For Q=b, the situation is somewhat more unstable, i.e. a given relative 
variation in sq is accompanied by roughly the same relative variation of Mp and fp- 



On the other hand, for a fixed value of s , the predicted M D and M B change by 15 % 
and 25 %, respectively, inside the sum rule windows in the Laplace parameter M 2 . In 
spite of all this, it is true that the final uncertainty in f P , due to the uncertainties in 
the input parameters, is smaller with the Laplace transform QCD sum rules than with 
the Hilbert moments. This can be appreciated from Tables 3 and 4, where I present the 
results obtained with the Laplace sum rules. However, the two types of sum rules exhibit 
different sensitivities to changes in the input parameters, in addition to having different 
advantages and shortcomings. For this reason they should be viewed as complementary 
methods within the general framework of QCD sum rules. 

Table 3: Laplace transform QCD sum rule results for Q = c 





m c = 


1.3 GeV 


m c = 


1.4 GeV 


Id/ fir 


1.30 


- 1.46 


1.15 


- 1.29 


s \d (GeV 2 ) 


5.5 


- 5.5 


5.0 


- 5.0 




1.56 


- 1.75 


1.45 


- 1.64 


s \d s (GeV 2 ) 


6.0 


- 6.0 


5.5 


- 5.5 


IdJ Id 


1.20 


- 1.20 


1.26 


- 1.27 



Table 4: Laplace transform QCD sum rule results for Q = b 





m b = 4.67 GeV 


m b = 4.77 GeV 




1.32 - 1.40 


1.14 - 1.18 


s \b (GeV 2 ) 


35.0 - 35.5 


34.5 - 35.0 


JbJ fn 


1.58 - 1.68 


1.38 - 1.44 


so\b s (GeV 2 ) 


35.5 - 36.0 


35.0 - 35.5 


JbJ Ib 


1.20 - 1.20 


1.21 - 1.22 



Given the fact that results from the Laplace sum rules have less of a spread than those from 
the Hilbert moments, one may be tempted to consider the former as the best determination 
of the leptonic decay constants. However, the two methods are not indepenedent, as the 
various vacuum condensates enter both sum rules, albeit with different weight factors, and 
different signs in the case of d=5 and d=6. Short of performing a correlation analysis, I 
feel one should not discard the results from the Hilbert moments (nor perform any average 
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from the two methods), but rather read the absolute minimum and maximum values from 
both sum rules in conjunction. In any case, the predictions for the ratios /d 8 //d, and 
fs a /f B turn out to be far less dependent on the values of the heavy quark masses, and 
the particular sum rule, leading to the accurate and stable predictions 

f D ./f D = 1.21 ±0.06, (12) 

f B Jf B = 1-22 ±0.02. (13) 
These results are in nice agreement with the expectation that the ratio between (12) and 



(13) should be close to unity (see e.g. ]I2"|| ). 

A comparison of the results listed in Tables 1-4 with predictions from lattice QCD (see 
e.g. fll3f ) shows reasonable agreement for fp, and fn a , but not quite for f B , although the 
ratio f Bs / f B does compare well. In the framework of fully relativistic QCD sum rules, it 
is simply not possible to obtain values of f B bigger than what is shown in Tables 1 - 4, if 
one uses the current best values of the input parameters. There is a recent claim to the 



contrary by Narison |I4" |, but I have been unable to reproduce his results, which I believe 
to be incorrect. 

It is only when one considers the infinite quark mass limit, and after resumming the large 
logarithms [15j , that one can approach lattice QCD predictions. The price to pay, though, 



is a two - loop correction at the 100 % level in the expression for f\. 
Acknowledgements 

The author wishes to thank the organizers of the Third Workshop on the Tau-Charm 
Factory for their effort in making this conference a great success. He also thanks Daniele 
Amati for his kind hospitality at SISSA, where this report was written. This work was 
supported in part by the Foundation for Research Development. 



References 



7 



[1] T.M. Aliev, V.L. Eletsky, Sov. J. Nucl. Phys. 38 (1983) 936. 

[2] C.A. Dominguez, N. Paver, Phys. Lett. B197 (1987) 423; B199 (1987) 596 (E). 

[3] L.J. Reinders, Phys. Rev. D38 (1988) 423. 

[4] S. Narison, Phys. Lett. B198 (1987) 104. 

[5] P. Colangelo, G. Nardulli, A. A. Ovchinnikov, N. Paver, Phys. Lett. B269 (1991) 
201. 

[6] C. A. Dominguez, E. de Rafael, Ann. Phys. (N.Y.) 174 (1987) 372; C. A. Dominguez, 
C. van Gend, N. Paver, Phys. Lett. B253 (1991) 241. 

[7] C. A. Dominguez, N. Paver, Phys. Lett. B253 (1991) 241. 

[8] R. A. Bertlmann et al, Z. Phys. C39 (1988) 231; C. A. Dominguez, J. Sola, Z. Phys. 
C40 (1988) 63. 

[9] D.J. Broadhurst, S.C. Generalis, Open University Report No. OUT-4102-8/R (1982); 
S.C. Generalis, Ph.D. Thesis, Open University Report No. OUT-4102-13 (1984). See 
also K. Schilcher, M.D. Tran, N.F. Nasrallah, Nucl. Phys. B181 (1981) 91. 

[10] C. A. Dominguez, N. Paver, SISSA Report No. SISSA Ref. 133/93/EP (1993). 

[11] ALEPH Collaboration, D. Buskulic et al, Phys. Lett. B311 (1993) 425. 

[12] B. Grinstein, SSCL Report No. SSCL - 492 (1993). 

[13] G. Martinelli, these proceedings. 

[14] S. Narison, Z. Phys. C55 (1992) 671. 

[15] E. Bagan et al, Phys. Lett. B278 (1992) 457; M. Neubert, Phys. Rev. D45 (1992) 
2451. 



8 



